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ABSTRACT 

This is an addendum to the paper by Cappellari (2008, MNRAS, 390, 71), which presented 
a simple and efficient method to model the stellar kinematics of axisymmetric stellar sys- 
tems. The technique reproduces well the integral-field kinematics of real galaxies. It allows 
for orbital anisotropy (three-integral distribution function), multiple kinematic components, 
supermassive black holes and dark matter. The paper described the derivation of the projected 
second moments and we provided a reference software implementation. However only the 
line-of-sight component was given in the paper. For completeness we provide here all the six 
projected second moments, including radial velocities and proper motions. We present a test 
against realistic N-body galaxy simulations. 
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1 PROJECTED SECOND MOMENTS 

In lCappellari] d2008h we used the lJeanj il922h equations to derive 
the projected second velocity moments for an anisotropic (three- 
integral) axisymmetric stellar system w ith the density describ ed via 
the Multi-Gaussian Expansion (MGE. lEmsellem et al1ll994h . We 
stated in note 5 that the components of the proper motion disper- 
sion tensor can be can be written via single quadratures without the 
need for special functions, and provided a reference software im- 
plementation, called the Jeans Anisotropic MGE (JAM) methocQ. 
However only the line-of-sight component was given in the paper. 
For completeness we provide all expressions in this addendum. 

We a d opt id entical notation and coordinates system as in 
Cappellari (2008), and we refer the reader to that paper for de- 
tails and definitions. Following the approach out lined in note 5 of 
the pa per, we use the general formula^ (A5) of lEvans & de Zeeuwl 
( 1994) to write any of the six projected second moments as 



J ° J=lj=l 

exp{-A[x' 2 + y /2 (A + B)/£]} 



2 t 



du, 



(1-Cu*) j£ [l-(l-g|) W 2] 



(1) 



where a and (3 stand for any of the three projected coordinates x', 
y and z', and we defined 



£ = A + Bcos 2 
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2 With the substitution x, y — ¥ y,x 



(2) 



The line-of-sight ve locity second moment v- 



v z iv z i given 

in eq. (28) of Cappellari d2008h is obtained from equation (fTJi with 



T- 2 2/ 2 • , j - 3 A , rfi /3 ■ 2 ■ 

■F z'z' — a k1k ( cos i + OfcSin i)+T>x sin i. 



(3) 



The second moment v x iv x i of the proper motion in a direction par- 
allel to the projected galaxy major axis is obtained with 

= bkxrlql +V{[y' cos i (A + B)/£] 2 + sin 2 t/(2£)}. (4) 

The second moment v y iv y i of the proper motion in a direction par- 
allel to the the projected symmetry axis is obtained from T y i y i = 
J'z'z'i't /2 — i), which implies 

-,- 22/. 2-., 2 A . t> /2 2 ■ 

J- y i y i = Vklk ( sm I + Ok COS Ij+UX COS I. 



The expressions for the cross terms are 
Fx'y' = —T>x'y cos 2 i(A + B)/£, 



-T x i y i tanz = Dx y sini cosi (A + B)/£, 



T y i z ' = sini cosi [o k q k (1 — b/,) — T> x' 



(5) 

(6) 
(7) 
(8) 



The expressions for v x >v x ' and v y iv y i where also recently given in 
|P'Souza&Rb3 l l2012h . who presented an application of the JAM 
method to proper motions measurements. 



2 TESTING JAM AGAINST N-BODY SIMULATIONS 

To test the formalism of the previous section we compared the six 
second moments predicted by the JAM method against the ones 
measured by direct summation f rom the N-body parti cles of two 
realistic galaxy simulations from lLablanche et"aH d2012h . to which 
we refer the reader for details. Inside each projected pixels, with 
coordinates (x[, y' p ) and iV particles, the second moments are 
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Figure 1. Projected second velocity moments of simulations versus JAM models. Top Row: Voronoi-binne d and linearly in t erpola ted second moment 
calculated by direct summation (equation (5J) of the N-body particles of the evolved simulations N4179axi of Lablanche et al. (2012). The simulation is 
projected at an inclination i = 60° . The corresponding velocity moment is written at the top of each image. Second Row: Same as in the top row for the second 
moments predicted by the JAM model (equation {T}) at the known inclination i = 60° . The MGE fit to the surface brightness of the above simulation was 
used as input to the JAM procedure. The anisotropy /3 Z was fixed to the known average value /3f IM from the simulations. The model is a prediction and not 
a fit to the data: it has no free parameters! Third Row: Same as in the top row, for the simulation N4570axi, projected at i = 30°. Bottom Row: Same as in 
the second row for the JAM model based on the MGE fitted to the simulation of N4570axi, projected at i = 30°. The simulations do not precisely satisfy the 
JAM assumptions, but the simple JAM models, with no free parameters, are able to provide a striking prediction for all six second moments, for both realistic 
N-body simulations. 



V*V0(x' l ,y' p ) = —^(VajVpj). (9) 
3=1 

To reduce th e noise on the moments ma ps, we used the Voronoi 
binning method ( Cappellar i"&Copiil 2003) to group pixels in such 
a way that every bin contains about 500 N-body particles. We then 
used the software of Cappellari (2002) to fit an MGE model to the 
surface brightness of the simulation. The MGE parameters were 
used as input for the JAM model (equation (|T}) to predict all six 
projected second moments. The model is not fitted to the data and 
it has no free parameters! We fixed the inclination i to the known 
value and the anisotropy B* to the average value s from the simula- 
tion particles (Table CI of lLablanche et alj|2012h . 

Considering the simplicity of the model assumptions and the 
lack of free parameters in our test, the agreement between the re- 
alistic simulations and the JAM model is striking, confirming the 
usefulness of the approach. This suggest that the method may be 



applicable to the interpretation of the observations from future gen- 
erations of proper motions surveys (generally using a maximum 
likelihood approach for discrete data). It can be especially useful to 
break the degeneracies that affect some of the quantities obtained 
from projected kinematics alone (e.g. the dark matter profiles). 
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